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We  present  density-functional  theory  calculations  of  the  chemisorption  of  atomic  species  0,  S,  C,  H  and 
reaction  intermediates  OH,  SH,  and  CHn  (n  =  1,  2,  and  3)  on  M/Ni  alloy  model  catalysts  (M  =  Bi,  Mo,  Fe,  Co, 
and  Cu).  The  activity  of  the  Ni  alloy  catalysts  for  solid-oxide  fuel  cell  (SOFC)  anode  oxidation  reactions  is 
predicted,  based  on  a  simple  descriptor,  i.e.,  the  binding  energy  of  oxygen.  First,  we  find  that  the  binding 
of  undesirable  intermediates,  such  as  C  and  S,  can  be  inhibited  and  the  catalytic  activity  of  planar  Ni-based 
anodes  can  be  tuned  towards  oxidation  by  selectively  forming  a  bimetallic  surface  alloy.  In  particular, 
Cu/Ni,  Fe/Ni,  and  Co/Ni  anode  catalysts  are  found  to  be  most  active  towards  anode  oxidation.  On  the 
other  hand,  the  Mo/Ni  alloy  surface  is  predicted  to  be  the  most  effective  catalyst  in  terms  of  inhibiting 
the  deposition  of  C  and  S  (while  still  preserving  relatively  high  catalytic  activity).  The  formation  of  a 
surface  alloy,  which  has  the  alloy  element  enriched  on  the  topmost  surface,  was  found  to  be  critical  to 
the  activity  of  the  Ni  alloy  catalysts. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid-oxide  fuel  cells  (SOFCs)  have  emerged  as  a  promising  tech¬ 
nology  for  efficient  electrochemical  energy  conversion.  In  addition, 
their  fuel  flexibility  enables  the  direct  utilization  of  natural  gas 
by  internal  steam  reforming  [1,2].  A  typical  SOFC  contains  four 
basic  components:  anode,  electrolyte,  cathode,  and  interconnect. 
For  many  years,  yttria-stabilized  zirconia  (YSZ)  has  been  commonly 
used  as  the  electrolyte,  and  a  Ni/YSZ  cermet  has  served  as  the 
anode,  which  is  where  the  fuel  is  electrocatalytically  oxidized  at 
the  three-phase  boundary  (TPB)  of  the  fuel/electrode/electrolyte. 
SOFCs  typically  operate  at  high  temperatures  (in  the  range  of 
800-1 000  °C),  which  enables  high  ionic  conductivity  (O2-  diffu- 
sivity)  within  the  YSZ  lattice,  superior  catalytic  activity  of  the 
Ni-based  anode  towards  oxidation  of  the  fuel  (mainly  H2  and  CO), 
and  steam  reforming  of  hydrocarbons.  However,  the  high  operating 
temperatures  also  cause  severe  materials  degradation,  application 
limitations,  and  high  integration  costs. 

The  catalytic  performance  of  the  anode  in  a  SOFC  plays  a  key 
role  in  determining  the  overall  fuel  cell  performance.  The  devel¬ 
opment  of  optimal  anode  materials  (i.e.,  combinations  of  electrode 
and  electrolyte)  with  high  catalytic  activity  at  lower  temperatures 


*  Corresponding  author.  Tel.:  +1  205  348  1733;  fax:  +1  205  348  7558. 
E-mail  address:  hturner@eng.ua.edu  (C.H.  Turner). 


(i.e.,  500-800  °C)  has  been  an  active  research  topic  [3-11  ].  Although 
a  number  of  material  candidates  have  been  considered  (see  lat¬ 
est  review  [12]  and  references  therein),  Ni-cermet-based  materials 
still  remain  the  most  efficient  and  practical  anode  materials.  This 
is  in  large  part  due  to  the  fact  that  the  basic  properties  of  such 
‘standard’  Ni-zirconia  or  Ni-ceria  cermets  are  well-established.  In 
particular,  they  possess  superior  catalytic  activity,  electrical  con¬ 
ductivity,  and  gaseous  diffusivity,  as  well  as  excellent  structural 
integrity  under  operating  conditions.  Hence,  further  development 
of  Ni-cermet-based  anodes  should  be  realized  by  fabricating  opti¬ 
mal  microstructures  and  by  identifying  moderate  compositional 
optimization  of  existing  Ni-cermet  systems  [3-9].  By  pursuing  this 
route,  the  inherent  catalytic  functionality  and  compatibility  with 
the  other  components  in  the  fuel  cell  and  the  fuel  cell  stack  should 
be  conserved. 

One  of  the  main  advantages  of  Ni-cermet-based  anodes  is  the 
high  electrocatalytic  activity  of  Ni  towards  H2,  CO,  hydrocarbon 
oxidation,  and  catalytic  activity  of  Ni  towards  steam  reforming  of 
methane  (SRM)  under  operating  conditions.  However,  Ni  also  cat¬ 
alyzes  the  formation  of  carbon  from  hydrocarbons  under  reducing 
conditions  [13],  which  can  block  active  Ni  sites.  In  addition,  the  Ni 
catalyst  is  easily  poisoned  by  sulfur  (a  common  impurity  in  natu¬ 
ral  gas),  which  can  also  block  active  Ni  sites  and  possibly  diffuse 
inside  nickel.  In  certain  cases,  Ni-S  eutectic  liquids  can  even  be 
formed  under  operating  conditions  [14].  Ultimately,  carbon  depo¬ 
sition  and  sulfur  poisoning  will  lead  to  the  deactivation  of  the  anode 
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Fig.  1.  Optimized  (111)  surfaces  with  various  3-fold  hollow  sites  labeled:  (a)  Ni,  (b)  Nio.75Cuo.25,  and  (c)  Nio.sCuo.s.  Blue  =  Ni  and  orange  =  Cu.  (For  interpretation  of  the 
references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


and  loss  of  its  durability.  Therefore,  further  development  requires 
Ni-based  anodes  with  enhanced  tolerance  to  coking  and  sulfur  poi¬ 
soning.  With  respect  to  the  selection  of  alternate  anode  materials, 
many  other  transition  metals  besides  Ni  have  been  investigated; 
however,  none  of  them  are  able  to  compete  with  and  replace  Ni 
completely,  in  terms  of  performance  and  cost  [12,15,16].  On  the 
other  hand,  Ni-alloys  can  help  mitigate  carbon-  and  sulfur-induced 
deactivation,  while  maintaining  the  inherent  material  structure 
and  performance  of  the  SOFC  [17-22].  For  instance,  carbon  for¬ 
mation  can  be  suppressed  by  employing  Ni  surface  alloys,  such 
as  Cu/Ni(l  1 1)  [22,23]  or  Au/Ni(l  1 1)  [24].  Furthermore,  several 
experimental  studies  of  Ni-alloy  anodes  have  demonstrated  the 
beneficial  effects  of  the  alloying  element  (e.g.,  Mo,  Fe,  Cu,  Pd,  and 
Sn)  on  the  fuel  cell  performance  [13,17-21]. 

In  order  to  provide  additional  insights  into  the  development  of 
Ni-alloy  anode  materials,  we  employed  density  functional  theory 
(DFT)  to  study  the  effects  of  different  alloying  elements  on  the  cat¬ 
alytic  properties.  Mainly,  we  predict  the  binding  energies  of  various 
adsorbates  and  reaction  intermediates  of  SOFC  anode  reactions  on 
bimetallic  M/Ni(l  1 1 )  surface  alloys  (M  =  Bi,  Mo,  Fe,  Co,  and  Cu).  We 
calculate  the  adsorption  of  atomic  species  O,  S,  C,  and  H  and  reac¬ 
tion  intermediates  OH,  SH,  and  CHn  (n  =  1,  2,  and  3)  on  Ni  surface 
alloys,  with  different  alloying  compositions,  and  this  information 
is  used  to  quantify  the  expected  catalytic  activity  for  oxidation  at 
the  anode.  We  find  that  the  binding  of  undesirable  intermediates 
(C  and  S)  can  be  suppressed,  and  the  catalytic  activity  of  the  Ni- 
alloy  anodes  can  be  tuned  towards  oxidation  by  forming  selective 
bimetallic  M/Ni(l  1 1)  surface  alloys. 


2.  Methods 

All  calculations  were  performed  using  DFT,  as  implemented 
in  the  Vienna  ab  initio  simulation  package  (VASP)  [25].  The  pro¬ 
gram  uses  a  plane-wave  basis  set,  thus  offering  good  access  to  the 
Hellmann-Feynman  forces  acting  on  all  of  the  atoms  in  the  super¬ 
cell.  The  spin-polarized  GGA  PBE  functional  [26],  an  all-electron 
plane-wave  basis  set  with  an  energy  cutoff  of  400  eV,  and  the 
projector  augmented  wave  (PAW)  method  [27,28]  were  adopted. 
The  Brillouin-zone  of  the  p( 2  x  2)  lateral  supercell  was  sampled 
by  5  x  5  x  1  k-points  using  the  Monkhorst-Pack  scheme,  and  first- 
order  Methfessel-Paxton  smearing  [29]  of  0.2  eV  was  employed  in 
the  integration  to  speed  up  the  convergence.  The  optimization  was 
performed  using  conjugate-gradient  algorithm  without  symmetry 
constraints,  and  the  convergence  threshold  was  set  to  be  10-4  eV 
for  total  energy  and  10-2  eVA-1  in  force  on  each  atom.  All  of  the 
reported  energies  were  extrapolated  to  kBT  =  0  eV. 

The  Ni(  111)  surface  was  modeled  by  a  four-layer  slab,  with  the 
lower  two  layers  fixed  at  their  equilibrium  bulk  phase  positions 


(calculated  lattice  constant  =  3.522  A),  while  the  upper  two  layers 
were  allowed  to  relax.  The  two  successive  slabs  were  separated  by 
a  13  A  vacuum  region  to  ensure  that  the  adsorbates  (i.e.,  O,  S,  C, 
H,  OH,  SH,  and  CHn,  n  =  1 ,  2  and  3)  and  the  subsequent  slab  would 
not  interact.  The  M/Ni(l  1 1)  surface  alloy  (M  =  Bi,  Mo,  Fe,  Co,  and 
Cu)  was  modeled  in  the  same  fashion  as  the  pure  Ni(l  11)  surface. 
However,  in  the  alloy,  the  Ni  atoms  were  substituted  with  M  atoms 
on  the  topmost  layer  in  the  optimized  p( 2  x  2)  lateral  supercell, 
with  four  metal  atoms  in  each  of  the  four  layers  of  the  slab  supercell, 
as  shown  in  Fig.  1.  The  M/Ni(l  1 1 )  surface  alloys  were  constructed, 
with  6=  1/4  and  1/2  ML  (in  terms  of  M  coverage)  by  substituting 
one  or  two  of  the  surface  nickel  atoms  with  M  atoms.  Alternatively 
for  copper,  Cu/Ni(l  1 1 )  surfaces  were  also  cleaved  from  a  NiCu  bulk 
alloy  having  one,  Fig.  1(b),  or  two,  Fig.  1(c),  copper  atoms  in  each 
layer. 

The  binding  energy  (BE)  in  this  work  is  defined  as 
BE  =  Etotal  [  slab ]  +  Etotal  [ adsorbate ]  -  Etotal  [  slab  +  adsorbate ] ,  where 
Etotai  is  the  total  energy  of  the  system  in  one  supercell,  and  the 
gas-phase  adsorbate  is  used  as  the  reference  for  £total[ adsorbate]. 

3.  Results  and  discussion 

In  general,  the  net  SOFC  electrode  reactions  can  be  written  as 
(all  molecules  are  in  the  gas  phase  and  all  atomic  and  intermediate 
species  are  in  adsorbed  states): 

Anode  :  H2 +  02_  (electrolyte)  ->  H20  +  2e_  (1) 

CO  +  O2-  (electrolyte)  ->  C02  +  2e-  (2) 

Cathode:  02  +4e- -*  202- (electrolyte)  (3) 

Overall  reaction  :  CO  +  H2  +  02  -*  C02  +H20  (4) 

If  natural  gas  is  used  as  the  fuel,  the  internal  steam  reforming  of 
hydrocarbons  (mainly  methane)  also  occurs  at  the  anode: 

CH4  +  H20  -*  CO  +  3H2  (5) 

The  overall  reactions  in  (4)  and  (5)  consist  of  a  number  of  elemen¬ 
tary  reactions  involving  atomic  species  such  as  0,  S  (as  an  impurity 
in  natural  gas),  C,  and  H,  and  reaction  intermediates  such  as  OH, 
SH,  and  CHn  (n  =  1,  2  and  3).  Therefore,  the  binding  nature  of  these 
atomic  species  is  relevant  to  the  catalytic  performance  and  deac¬ 
tivation  of  the  SOFC  anode.  According  to  the  “Sabatier  principle”, 
a  good  catalyst  has  active  sites  which  bind  reaction  intermediates 
with  an  intermediate  strength.  In  other  words,  binding  that  is  too 
weak  or  too  strong  leads  to  a  catalyst  which  is  either  less  active 
or  easily  poisoned.  In  particular,  the  chemical  interaction  between 
the  atomic  O  and  the  anode  is  believed  to  play  a  critical  role  in 
the  overall  fuel  cell  performance  [15,16].  It  should  be  noted  that 
the  fundamental  mechanisms  associated  with  the  electrochemistry 
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Fig.  2.  Maximum  binding  energy  (BE)  of  atomic  species  on  Ni(1_X)Mx  (111)  surface  alloys,  as  a  function  of  mole  fraction  (x):  (a)  0,  (b)  S,  (c)  C  and  (d)  H.  Black,  blue,  purple, 
green,  and  red,  denote  M  =  Bi  (■),  Mo  (^  ),  Fe  (▼  ),  Co  (  ),  and  Cu  (A  and  A  ),  respectively.  Hollow  red  triangles  represent  the  surface  cleaved  from  a  stoichiometric 

Ni(!_X)Cux  bulk  alloy.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


in  reaction  (4)  are  still  not  completely  understood  [30].  Two  main 
reaction  mechanisms  at  the  anode  have  been  widely  proposed:  (i) 
oxygen  spillover;  (ii)  hydrogen  spillover.  In  the  former  pathway, 
O2-  migrates  from  the  YSZ  electrolyte  to  the  Ni  surface  forming 
adsorbed  O  species,  which  oxidize  the  fuel  (i.e.,  CO  and  H2): 

O2-  (electrolyte)  ->  O  +  2e_  (6) 

In  the  latter  pathway,  H2  diffuses  through  the  porous  Ni-cermet 
forming  dissociated  H  atoms  on  the  Ni  surface,  which  subsequently 
release  their  electrons  to  the  Ni  electrode  and  hop  to  the  O2-  or  OH- 
site  on  the  ceramic,  YSZ,  forming  water.  The  neutralization  (charge 
interactions)  of  the  acidic  H+  and  the  basic  O2-  or  OH-  on  YSZ  serves 
as  the  driving  force  for  proton  migration.  In  both  pathways,  the 
electrons  (without  passing  through  the  ion  conductive  electrolyte) 
are  directed  towards  the  cathode  through  an  external  circuit,  which 
creates  electrical  current.  The  “oxygen  spillover”  mechanism  also 
explains  the  electrochemical  oxidation  of  CO  without  the  involve¬ 
ment  of  H2.  Thus,  we  use  the  binding  of  atomic  O  on  the  metal  as 
a  key  descriptor  for  the  performance  of  the  proposed  M/Ni  alloy 
catalysts. 

3.1.  Adsorption  of  atomic  species 

For  the  anode  oxidation  to  proceed,  the  fuel,  i.e.,  hydrogen  and 
hydrocarbons,  must  be  chemisorbed  and  dissociated  into  atomic 
species,  which  must  further  react  with  adsorbed  0  at  the  TPB.  We 
first  investigated  the  adsorption  of  atomic  species  0,  S,  C,  and  H  on 
various  M/Ni(l  1 1)  alloys  (M  =  Bi,  Mo,  Fe,  Co,  and  Cu).  We  began 
these  calculations  by  first  investigating  the  influence  of  the  surface 
composition  on  the  adsorption,  versus  the  influence  of  the  bulk 
catalyst  composition. 

Fig.  2  shows  the  maximum  binding  energy  (BE)  of  these 
atomic  species  as  a  function  of  the  mole  fraction  (x)  of  M  on 
the  Ni(1_x)Mx(l  1 1)  alloy  surfaces.  We  find  that  the  most  stable 
adsorption  sites  are:  (1)  the  [Ni— Ni— Ni]  3-fold  hollow  sites  on 
pure  Ni(l  1 1)  and  on  the  M/Ni(l  1 1)  alloy  surfaces,  with  0.25  ML 
coverage  of  M;  and  (2)  the  [Ni-Ni-M]  3-fold  hollow  sites  on  the 
M/Ni(l  1 1)  alloy  surface  with  0.5  ML  coverage  of  M  (see  Fig.  1). 
From  Fig.  2,  one  can  see  that  the  incorporation  of  Bi  has  the  most 
significant  effect  on  reducing  the  binding  strength  of  the  O,  S,  C,  and 
H  atomic  adsorbates.  However,  according  to  our  calculated  results, 
Bi  is  too  large  [31]  and  segregates  above  the  surface  if  two  Ni  and 
two  Bi  atoms  are  in  the  surface  layer.  This  only  allows  a  maximum 
Bi  coverage  of  0.25  ML  on  a  Bi/Ni(l  1 1)  surface,  in  our  modeling 
scheme,  where  each  surface  Ni  atom  has  four  surface  Ni  atoms  and 
two  surface  Bi  atoms  as  neighbors.  In  comparison,  Mo,  Fe,  Co,  and 


Cu  (being  smaller)  can  all  form  stable  surface  alloys  at  higher  con¬ 
centrations.  At  the  higher  0.5  ML  coverage  of  M,  each  Ni  atom  has 
two  Ni  atoms  and  four  M  atoms  as  neighbors,  and  equivalently  each 
M  atom  has  four  Ni  atoms  and  two  M  atoms  as  neighbors. 

Comparatively,  the  Mo/Ni(l  1 1 )  surface  alloy  with  0.25  ML  cov¬ 
erage  of  Mo  has  the  best  performance  in  terms  of  lowering  the 
binding  strength  of  C  and  S,  while  maintaining  the  binding  strength 
of  0  at  a  favorable  level,  which  ensures  a  high  catalytic  activity 
towards  oxidation  (see  Section  3.3).  Our  calculated  results  agree 
very  well  with  the  early  experimental  finding  that  doping  the  Ni 
anode  with  small  quantities  of  molybdenum  (<1%)  leads  to  a  sig¬ 
nificant  reduction  in  carbon  deposits,  while  having  little  effect  on 
the  SRM  and  SOFC  performance  [13].  If  the  Mo  coverage  is  fur¬ 
ther  increased  up  to  0.5  ML,  there  is  a  significant  increase  in  the 
binding  strength  of  O  and  a  slight  increase  in  the  binding  strength 
of  S.  Both  values  are  higher  than  on  a  pure  Ni(  111)  surface.  The 
beneficial  effect  of  Mo  on  inhibiting  C  adsorption  remains  intact, 
regardless  of  the  higher  Mo  surface  coverage.  By  contrast,  the  bind¬ 
ing  strength  of  H  increases  on  the  Ni  alloy  with  0.25  ML  Mo,  but  the 
binding  weakens  at  a  coverage  of  0.5  ML  of  Mo,  as  compared  to  that 
on  a  pure  Ni(l  1 1)  surface.  Since  the  “hydrogen  spillover”  mech¬ 
anism  directly  depends  on  the  chemisorption  of  H  atoms,  which 
are  further  split  into  protons  and  electrons,  this  alternative  reac¬ 
tion  pathway  can  be  sensitive  to  the  binding  of  H  atoms  on  Ni  the 
alloy. 

Likewise,  the  Fe/Ni  and  Co/Ni  alloys  demonstrate  behavior  sim¬ 
ilar  to  the  Mo/Ni  alloy,  in  terms  of  modifying  the  binding  energy, 
but  with  a  more  moderate  influence.  The  trend  predicted  from  the 
Cu/Ni  alloy  surface  is  different  from  that  of  Bi,  Mo,  Fe,  and  Co.  From 
Fig.  2,  one  can  see  that  the  binding  strengths  of  0,  S,  C,  and  H  are 
all  enhanced  at  0.25  ML  of  Cu  and  then  weaken  at  0.5  ML  of  Cu,  as 
compared  to  those  on  a  pure  Ni(  111)  surface.  This  suggests  that 
the  beneficial  effect  of  Cu  on  inhibiting  S  and  C  is  more  pronounced 
at  higher  Cu  coverage,  which  agrees  well  with  the  experimental 
observations  that  Cu/Ni  catalysts,  with  Cu  concentrations  up  to 
50%,  have  demonstrated  the  highest  ratio  of  reforming:  coking  rates 
[32,33  ].  Our  previous  study  also  showed  that  Cu  addition  can  inhibit 
carbon  formation  during  the  sequential  dissociation  of  CH4,  as  the 
activation  energy  barrier  of  the  final  step  is  increased  [23]. 

As  a  comparison,  we  also  investigated  the  adsorption  behavior 
ofO,  S,  C,  and  Hon  Cu/Ni(l  1 1)  alloy  surfaces  [see  Fig.  1(a)  and  (b)], 
cleaved  from  a  Cu/Ni  bulk  alloy  yielding  the  dashed  lines  in  Fig.  2. 
One  can  see  that  the  binding  strength  trends  of  the  atomic  species 
0,  S,  C,  and  H  on  the  bulk  Cu/Ni(l  1 1 )  alloy  surface  are  very  similar 
to  the  values  calculated  on  the  Cu/Ni(l  1 1)  surface  alloy  (i.e.,  the 
model  with  the  alloy  element  only  present  in  the  top  layer).  This 
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Fig.  3.  Binding  energy  (BE)  of  fragments  on  (1 1 1)  surfaces  as  a  function  of  the 
number  (z)  of  Cu  atoms  at  hep  (black)  and  fee  (green)  3-fold  hollow  sites,  as  shown 
in  Fig.  2,  wherez  =  0, 1, 2  corresponds  to  (1 1 1)  surfaces  of  Ni,  Nio.75Cuo.25,  Nio.sCuo.s 
bulk  alloy,  respectively.  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  the  article.) 


suggests  that  the  alloy  atoms  in  the  topmost  layer  play  the  largest 
role  in  mediating  the  adsorbate  binding  strength,  and  thus  the  cat¬ 
alytic  properties  of  the  surface.  Therefore,  it  is  critical  to  know  the 
precise  surface  alloy  composition  in  the  experimental  systems.  Fur¬ 
thermore,  this  provides  supporting  evidence  that  precious  metal 
loading  can  be  minimized,  if  the  surface  alloy  composition  can  be 
maintained  at  the  desired  level. 

As  shown  in  Fig.  3,  the  BEs  of  0,  S,  C,  and  H  decrease  almost 
linearly  as  the  number  (z)  of  Cu  atoms  around  the  3-fold  hollow 
sites  increases.  This  trend  continues,  regardless  of  the  particular 
hep  or  fee  site  configuration,  where  z  =  0, 1 ,  and  2  corresponds  to  Ni, 
Nio.75Cuo.25,  and  Ni0.5Cu0.5(l  1 1)  surfaces,  respectively.  This  illus¬ 
trates  the  important  effect  of  the  alloying  element  in  the  ensemble 
of  nearest-neighbor  atoms.  The  deposition  of  C  and  S  generally 
requires  large  Ni  ensembles.  By  contrast,  the  adsorption  of  O,  S, 
C,  and  H  is  enhanced  at  the  [Ni-Ni-Ni]  3-fold  hollow  sites  by  Cu 
alloying  (see  Fig.  2),  as  a  result  of  the  Cu  ligand  effect. 


3.2.  Adsorption  of  reaction  intermediates 

In  a  high  temperature  SOFC,  internal  steam  reforming  of 
methane  occurs  either  when  wet  natural  gas  is  imported  as  a  feed¬ 
stock  or  as  a  result  of  the  reaction  of  CFI4  with  the  FI20  produced 
from  H2  oxidation  at  the  anode.  It  has  been  suggested  that  ‘oxygen 
spillover’  is  the  dominant  reaction  mechanism  at  the  anode  [16].  In 
this  process,  oxygen  ions  are  transferred  by  surface  spillover  from 
the  YSZ  electrolyte  to  the  Ni  surface,  followed  by  reaction  with  H 
and  CHn  (produced  from  the  dehydrogenation  of  CFI4).  Adsorbed 
atomic  O  species,  which  can  also  be  produced  from  the  dissocia¬ 
tion  of  FI20  on  the  Ni  surface,  play  a  key  role  in  the  SOFC  anode 
reactions.  At  high  operating  temperatures  (typical  of  SOFCs),  CFI4 
is  most  likely  to  be  completely  dehydrogenated  to  C,  unless  a  much 
higher  activation  energy  barrier  exists  to  inhibit  the  final  dehydro¬ 
genation  step:  CH  ->  C  +  H.  An  alternate  (and  preferred)  route  is  the 
oxidation  of  the  CFI  species:  CFI  +  0  ->  FI  +  CO  [31  ]. 

As  shown  in  Fig.  4,  we  predict  the  adsorption  of  other  reac¬ 
tion  intermediates  on  the  M/Ni(l  1 1 )  alloy  surfaces,  including:  CHn 
(n  =  1 , 2  and  3 ),  OH,  and  SH  fragments.  The  predictions  are  based  on 
the  calculated  binding  energy  of  A  and  using  a  scaling  relation¬ 
ship  to  extrapolate  the  BEs  of  AHn  (A  =  C,  O  and  S).  This  approach  is 
described  in  more  detail  in  Ref.  [32],  and  it  is  expected  to  be  appli¬ 
cable  to  Ni-based  surface  alloys  [23].  From  Fig.  4,  one  can  see  that 
when  OH  and  SH  are  adsorbed  at  the  same  sites  on  the  Ni(  111)  sur¬ 
face,  the  predicted  BEs  are  quite  close  in  energy.  The  same  behavior 
is  seen  on  the  M/Ni(l  1 1 )  alloy  surfaces,  suggesting  that  OH  and  SH 
are  competitive  reaction  intermediates  during  the  elementary  reac¬ 
tion  processes.  The  BEs  of  CHn,  OHn,  and  SHn  increase  as  n  decreases 
(for  C,  O  and  S  compare  with  Fig.  2),  and  this  trend  is  consistent  with 
the  increase  in  the  free  bonding  capability  of  the  AHn  species.  This 
increasing  BE  also  serves  as  the  driving  force  for  the  dissociation 
of  CH4,  H20,  and  H2S  in  the  fuel,  leading  to  the  formation  of  the 
atomic  species  (i.e.,  C,  O,  and  S)  on  the  catalyst  surface. 

3.3.  Catalytic  activity  of  bimetallic  M/Ni  surface  alloys  for  SOFC 
anode  reactions 

Previous  studies  have  shown  that  oxygen  spillover,  where 
adsorbed  O  is  a  key  intermediate,  can  be  a  dominant  reaction  chan¬ 
nel  under  typical  SOFC  operating  conditions  [15,16].  This  suggests 
that  the  BE  of  O  may  be  used  as  a  simple  descriptor  of  the  cat¬ 
alytic  activity  of  the  anode.  In  Fig.  5,  we  draw  the  characteristic 
volcano-type  [33]  plot  of  the  calculated  catalytic  activity  for  a  col¬ 
lection  of  anode  oxidation  reactions  versus  the  0  binding  strength 
on  (1 1 1)  transition  metal  surfaces,  with  lines  re-plotted  from  the 
study  by  Rossmeisl  and  Bessler  [16].  In  their  study,  Ni  was  iden¬ 
tified  as  the  best  catalyst  for  the  SOFC  anode  reactions  among  a 
number  of  transition  metals,  including  Ru,  Rh,  Pd,  Pt,  and  Au.  It  can 
be  clearly  seen  that  the  bimetallic  Ni-alloy  systems  studied  in  this 


Fig.  4.  Maximum  binding  energy  (BE)  of  CHn  (n  =  1,  2  and  3),  OH,  and  SH  fragments  on  Ni(i_x)Mx  (1 1 1)  surface  alloys,  as  a  function  of  the  mole  fraction  (x):  (a)  Bi,  (b)  Mo,  (c) 
Fe,  (d)  Co,  and  (e)  Cu. 
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Fig.  5.  The  volcano  curve  of  theoretical  activity  versus  binding  energy  (BE)  of  O  on 
(111)  transition  metal  surfaces,  re-plotted  from  Fig.  6  of  Ref.  [15]. 


work  possess  optimal  activity,  close  to  that  of  pure  Ni.  In  particu¬ 
lar,  Nio.sCuo.s,  Nio.75Coo.25.  Nio.75Feo.25.  Nio.sCoo.s,  and  Nio.75Cuo.25 
surface  alloys  have  superior  catalytic  activity  towards  anode  oxi¬ 
dation.  Thus,  improvements  can  be  made,  without  abandoning  the 
advantages  of  Ni  (low  cost,  material  compatibility,  etc.).  For  exam¬ 
ple,  our  predications  on  NiFe  surface  alloys  agree  well  with  a  very 
recent  study  [34],  which  shows  that  a  NiFe  bimetallic  alloy  used  as 
an  anode  can  significantly  improve  the  maximum  power  density 
and  lower  the  anodic  overpotential  at  intermediate  temperatures. 
Moreover,  the  addition  of  a  small  amount  of  Fe  to  Ni  can  drastically 
decrease  the  aggregation  of  Ni  under  the  cell  operating  conditions. 

4.  Conclusions 

The  ideal  materials  choice  for  SOFC  anodes  is  determined  by  a 
number  of  factors,  such  as  catalytic  activity,  electrical  conductivity, 
thermal  compatibility,  and  chemical  stability.  Thus,  the  forma¬ 
tion  of  Ni-alloy  surfaces  may  provide  advantages  that  improve  the 
performance  of  existing  Ni-cermet  based  anode  materials.  These 
efforts  may  lead  to  higher  efficiency  and  a  longer  lifetime  at  lower 
operational  temperatures.  In  this  work,  we  report  a  DFT-based 
computational  study  on  the  effects  of  alloying  elements  on  the 
adsorption  of  reaction  intermediates,  which  are  critical  for  quan¬ 
tifying  the  catalytic  activity  of  bimetallic  M/Ni(l  1 1)  surface  alloys 
(M  =  Bi,  Mo,  Fe,  Co,  and  Cu)  for  SOFC  anode  reactions.  We  find  that 
the  binding  of  undesirable  intermediates  (i.e.,  C  and  S)  can  be  sup¬ 
pressed  and  the  catalytic  activity  of  Ni-based  anodes  can  be  tuned 
towards  oxidation  by  selectively  forming  bimetallic  M/Ni(l  1 1) 
surface  alloys.  The  Cu/Ni,  Fe/Ni  and  Co/Ni  systems  appear  to  be  most 
active  towards  anode  oxidation,  based  on  the  correlation  between 
the  catalytic  activity  and  the  binding  energy  of  O.  On  the  other  hand, 
the  Mo/Ni  alloy  is  the  most  effective  catalyst,  in  terms  of  inhibiting 
the  deposition  of  C  and  S,  while  still  maintaining  a  high  catalytic 


activity.  The  formation  of  a  surface  alloy  is  found  to  be  critical  to 
tuning  the  activity,  and  this  is  particularly  important  if  precious 
metals  are  used  as  alloying  elements,  due  to  cost. 
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